The anisotropic Higgs oscillator on the two-dimensional sphere 
and the hyperbolic plane 



Angel Ballesteros 1 , Francisco J. Herranz 1 and Fabio Musso 2 

1 Departamento de Fisica, Universidad de Burgos, E-09001 Burgos, Spain 

2 Dipartimento di Fisica "Edoardo Amaldi", Universita Roma Tre, 1-00146 Rome, Italy 

E-mail: angelb@ubu.es, fjherranz@ubu.es, fmusso@ubu.es 



Abstract 

An integrable generalization on the two-dimensional sphere S 2 and the hyperbolic plane 
H 2 of the Euclidean anisotropic oscillator Hamiltonian with "centrifugal" terms given by 

n= 1 -( P l+pl) + s q l + (s + n) q 2 2 + ^ + ^ 

A HI '12 

is presented. The resulting generalized Hamiltonian H K depends explicitly on the Gaus- 
sian curvature k of the underlying space, in such a way that all the results here presented 
hold simultaneously for S 2 (k > 0), H 2 (k < 0) and E 2 (k = 0). Moreover, H K is explic- 
itly shown to be integrable (albeit not superintegrable) for any values of the parameters 
5, Q, Ai and A 2 . Therefore, H K can also be interpreted as an anisotropic generalization 
of the curved Higgs oscillator, that is recovered as the isotropic limit f2 — > of H K . Fur- 
thermore, numerical integration of some of the trajectories for ~K K are worked out and the 
dynamical features arising from the introduction of a curved background are highlighted. 
In particular we focus on the case Cl = 35, whose Euclidean limit n — > is the superin- 
tegrable 1:2 oscillator. In this case we illustrate how the Euclidean superintegrability of 
% is broken when the curvature k is non-zero. Since for the specific f2 = 35 case another 
superintegrable curved generalization is already known in the literature, the existence of 
an ^-dependent plurality of integrable curved systems whose Euclidean limit is the same 
anisotropic oscillator H can be conjectured. Finally, the geometric interpretation of the 
curved "centrifugal" terms appearing in H K is also discussed in detail. 
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1 Introduction 



The two-dimensional (2D) anisotropic oscillator with "centrifugal" (or "Rosochatius" ) terms 

n = Wi +Pl) + *Zi + (<* + tyq 2 2 + ^ + A |, (1.1) 

with (5, and Aj being real parameters and {qi,Pj} = 5y = 1,2), is one of the most 
elementary (albeit instructive) examples of completely integrable Hamiltonian systems. Its 
complete integrability (in the Liouville sense) is a straightforward consequence of its manifest 
separability in Cartesian coordinates, and for some specific choices of the parameters the 
system turns out to be superintegrable due to the existence of an additional independent 
integral of the motion. 

The aim of this paper is to present an integrable generalization of this Hamiltonian on 
the 2D sphere S 2 and on the hyperbolic plane H 2 , in the sense that the Hamiltonian (1.1) 
will be recovered when the Euclidean limit is considered. In fact, throughout the paper we 
shall make use of the constant Gaussian curvature k of these two curved spaces as an explicit 
parameter, and all the results here presented will be simultaneously valid for S 2 (k > 0), 
H 2 (k < 0) and E 2 (k = 0). In this way, the "flat /linear" limit k — > in the spherical 
and hyperbolic systems will be always well defined leading to the corresponding Euclidean 
Hamiltonian (1.1) (obviously, the term linear refers to the equations of motion when both Aj 
vanish). Conversely, the non-trivial dynamical nature of the transition from the "flat/linear" 
systems defined on E 2 to the "curved/nonlinear" ones associated to S 2 and H 2 can be clearly 
appreciated through this approach. 

Since integrability will be our essential guiding principle when constructing the curved 
generalization of the system, let us briefly recall the (Liouville) integrability properties of 
(1.1) in terms of the S, fi and A« parameters. Obviously, H is always integrable since it is 
separable in Cartesian coordinates: 

n=l 1 +l 2 , ll = lpl + 6q 2 1 + ^, l 2 = \p 2 2 + (5 + tt)q 2 2 + ^. (1.2) 

Thus each pair (T-L,Ti) and (H,X 2 ) is formed by two functionally independent functions in 
involution with respect to the canonical Poisson bracket. However, there are two particular 
situations for which (1.1) is superintegrable: 

• The isotropic oscillator (f2 = 0) with Rosochatius terms (Aj arbitrary), which is fre- 
quently called the Smorodinsky-Winternitz (SW) system [1, 2, 3, 4]. In this case the 
additional integral of the motion is quadratic in the momenta. 

• The anisotropic oscillator (Q ^ 0) with Rosochatius terms (A« arbitrary) provided that 
the ratio of the associated frequencies is a rational number [5, 6]. In this family of sys- 
tems only the commensurate 1:2 oscillator (fi = 3(5) plus a single Rosochatius potential 
(say A2 = 0) presents an additional integral which is quadratic in the momenta (if both 
Aj ^ or the anisotropic oscillator has any other pair of commensurate frequencies 
the additional integral is of higher order in the momenta [5, 6]). The resulting 1:2 
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superintegrable Hamiltonian and its quadratic integrals of motion turn out to be 



n' = \{p\+pl)+5ql+Uql 
2 at 




l' 2 = -p 2 2 +A5ql 



(1.3) 



£' = (qiP2 ~ qWi)Pi + 25(7x92 




It can be straightforwardly checked that the sets (H', ![,£') and {%' ^T^C) are formed 
by three functionally independent functions. 

To the best of our knowledge, there are only two particular cases of the Euclidean oscillator 
% (1.1) whose integrable generalization to the 2D sphere and hyperbolic plane are known. 
The first one is a curved generalization of the isotropic O = oscillator, which is the so-called 
Higgs oscillator [7, 8]; in this case when both centrifugal terms are taken different from zero 
the resulting Hamiltonian provides the curved SW system [9, 10, 11, 12, 13]. The second one 
is a curved analogue %' K of the superintegrable nonisotropic 1:2 system %' (1.3), which was 
found in [11]. Remarkably enough, these two curved systems are also superintegrable. 

Therefore, a "generic" (for any value of Q) integrable generalization of H (1.1) on S 2 and 
H 2 was still lacking, and we present it in this paper. We will denote it H K and, since the limit 
$7 — > of Tl K will lead to the Higgs oscillator, we shall call this new system the "anisotropic 
Higgs oscillator" . 

We would like to emphasize that % K seems to be superintegrable only in the Q = case 
(i.e., the isotropic Higgs oscillator or, equivalently, the 1:1 case). In fact, numerical analysis 
of bounded trajectories for H K with Q = 35 shows that they are not periodic, in contrast with 
what it should be expected for any superintegrable system. Moreover we find that the curved 
analogue of %' , i.e. %' K , is no longer a particular case of % K with Q = 35 and A2 = 0, which 
proves that both Hamiltonians define two different curved systems. On the contrary, bounded 
trajectories of H' R are found to be periodic in agreement with its superintegrability nature. 
All these results suggest that the curved superintegrable analogues of different commensurate 
Euclidean oscillators could be given by different systems that would "collapse" to a common 
form in the Euclidean limit k — > 0. 

The paper is organized as follows. In the next section we briefly recall the basics on Bel- 
trami and geodesic polar variables that generalize the Cartesian and polar ones, respectively, 
to S 2 and H 2 . In section 3 we introduce the anisotropic Higgs oscillator % K and we give the 
explicit form of its integral of the motion for generic values of the parameters. We also present 
some numerical analysis of its trajectories for some commensurate values of its Euclidean fre- 
quencies, showing that superintegrability does not survive under the deformation introduced 
through the curvature. In section 4 we analyse the superintegrable curved oscillator given by 
~H' K , and the 1:2 Lissajous curves on S 2 and H 2 are found as the trajectories for the system. 
In section 5 we perform a detailed geometrical analysis of the "curved" centrifugal terms [14] 
of all the previous curved oscillators, and their dynamical effects are numerically explored. 
Finally, the last section is devoted to some concluding remarks and open problems. 
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2 Geodesic polar and Beltrami variables 



First of all, let us rewrite the Euclidean Hamiltonians given in the previous section in terms of 
polar coordinates (r, <j>) with conjugate momenta (p r ,p<f,)- This change of coordinates induces 
the following canonical transformation 

, . , sin0 . cos0 
q 1 =rcos(p, q2 = rsm(p, pi=cos<pp r — p^, p 2 = sm(pp r -\ — P<t>, (2.1) 

that transforms the expressions (1.1)-(1.3) into 

H = I (pi + P \\ + 5r 2 + n r 2 sin 2 <f> + + X ' 



2 I r r 2 J " r 2 cos 2 <p r 2 sin 2 <j) ' 

1 / sin \ 2 A 

^•i = n ( cos0p r m ] +(5r 2 cos 2 0+^ , (2.2) 

2 \ r ) r l cos 2 

T 1 f • ^ , COS0 \ 2 2 . 2 A 2 

X 2 = - sm 4>p r H +(0 + ii)r sin + 



2 V r I r 2 sin <A 



^ = iU 2 + 4]+5r 2 + 3 ( 5r 2 sin 2 + 



2 \ rr r 2 / r 2 cos 2 0' 

1 / cos \ 2 
l{ = Xi , X 2 = - ( sin 4>p r H — + 45r 2 sin 2 0, 

„, / , sin0 \ o o 2Aitani 

£ = cos cppr Ps\P<i> + 2or cos q> sm - 

\ r / r cos <p 



(2-3) 



Now, in order to generalize the two above flat systems (2.2) and (2.3) to the sphere and to 
the hyperbolic plane we will apply a geometric aproach in which S 2 and H 2 are obtained as 
homogeneous spaces of certain 3D Lie groups. In this way, the realization of the associated 
symmetry generators in terms of geodesic polar and Beltrami (projective) coordinates will 
enable us to define the curved counterparts of these two systems. 



2.1 Geodesic polar coordinates 

Let so K (3) be the 3D real Lie algebra with generators Joi, J02 and J12 with commutation 
relations and Casimir operator given by 

[Jl2, Joi] = J02, [Jl2, J02] = — Joi, [Joi, J02] = KJl2, C = Jqi + Jq2 + K <^12) (2-4) 

where k is a real parameter. The three 2D classical Riemannian spaces of constant Gaussian 
curvature k can be collectively constructed as the homogenous spaces SO K (3)/ SO(2) where 
SO K (3) is the Lie group generated by so K (3) and SO{2) is the isotopy subgroup spanned by 
J\2- Hence Joi and J02 play the role of generators of translations, while J12 corresponds to 
the generator of rotations on the 2D space. In this way, according to the value of k we get 
the three particular homogeneous spaces: 

k > : Sphere k = : Euclidean plane k < : Hyperbolic plane 

S 2 = SO(3) /SO(2) E 2 = ISO(2) /SO{2) H 2 = SO{2, 1) /SO(2) 
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Figure 1 : Schematic representation of the symmetry generators Jij and geodesic distances of a point Q with 
geodesic polar coordinates (r, <f>) on a 2D homogeneous space. 



We recall that J12 leaves a point O invariant, the origin, and that each translation generator 
Joi moves O along a base geodesic k in such a manner that l\ and I2 are orthogonal at O as 
depicted in figure 1. 

These spaces can be globally embedded in the linear space M 3 with ambient or Weierstrass 
coordinates (xq,x) = (xo,xi,X2) by imposing the "sphere" constraint S: Xq + kx 2 = 1; the 
origin corresponds to the point O = (1,0) € M 3 . The geodesic polar coordinates r G (0, 00) 
and (f> 6 [0, 2tt) are defined as follows. If the particle is located at a point Q, then r is the 
distance between Q and the origin O measured along the geodesic / that joins both points. 
On the other hand, (j) is the angle which determines the orientation of / with respect to the 
base geodesic l\ (see figure 1). The ambient coordinates are parametrized in terms of (r,<f>) 
in the form [11, 15] 

xo = C K (r), xi = S K (r)cos0, x 2 = S K (r)sin0, 

where we have introduced the curvature-dependent functions 



C«(r) 



cos \fk r k > 
1 k = 

cosh J— k r k < 



S«(r) 



4= sin a/k ^ k > 

r k = 

—f= sinh V - « ^ k < 



The «;-tangent is defined by T K (r) = S K (r)/ C K (r). Properties for these ^-functions can be 
found in [11, 15, 16]. For instance, 



C£(r) + KS2(r) = l, 



Now let us denote by (p r , 
bracket for the three functions 



C K (r) = -fcS K (r), 



d 



S K (r) = C«(r) 



dr dr 
the conjugate momenta of (r, 0). The canonical Poisson 



sm< 



01 



COS <pp r 



T«(r) 



02 



COS (p 

sm0p r + - — 



T«(r) 



J12 = P<f> 



(2.5) 



gives us the Lie-Poisson analogue of the algebra so K (3) (2.4). Therefore, (2.5) provides a 
symplectic realization of the Lie-Poisson analogue of so K (3), and we can write the kinetic 
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energy T for the free motion of a particle moving on S 2 and H 2 as the above symplectic 
realization for the Casimir function of so K (3), namely: 



2 C " ^Ol + J 02 + Kj 12) - 2 1 -r ■ g 2 (r) 



T = = -(J 2 ! + J 2 2 + K jf 2 ) = -[p 2 r + ] • (2.6) 



Obviously, the n — > limit of this expression gives the 2D Euclidean kinetic energy. 
2.2 Beltrami coordinates 

On the other hand, if we apply the central projection with pole (0,0) G M 3 from (xo,x) to 
the Beltrami coordinates q£l 2 , that are defined as (0, 0) + fi (1, q) G S, we get the relations 

i q 

//=—====, x =n, x = /iq 



Vi + ^q 2 \A+ 



The image of this projection is the subset of M 3 with fi G R so that 1 + /-cq 2 > 0. In the 
hyperbolic or Lobachevsky plane with k = —\k\ < this gives the interior of the Poincare 
disk: 1 > |/-c|q 2 . Notice that when k = 0, the Beltrami coordinates reduce to the Cartesian 
ones on E 2 : xq = 1 and x = q. Hereafter we assume the following notation for q = (q±, q 2 ) 
and for the conjugate momenta p = (pi,P2)' 



q 2 = q\ + g| > p 2 = Pi + pi q • p = ?ipi + 92P2, |q| = \Jql + ql 

Then another symplectic realization of the Lie-Poisson algebra so K (3) is found to be [17] 

Joi =Pi + «(q- p)qi, i = l,2; hi = q\Pi - qiP\, (2.7) 

and the free Hamiltonian T that determines the geodesic motion on the Riemannian spaces 
is again obtained from the Casimir (2.4) under the above realization, namely 

T = 2° = 2 {J2qi + J2 ° 2 + ^ = 2 (1 + ^ ^ + K(q " p) *) • (2 - 8) 

We remark that, as expected, there is canonical equivalence between the Beltrami and 
geodesic polar variables through the canonical transformation given by 



2 sin (j) 



qi = T K (r)cos0, pi = C K (r) cos^p. 



T«(r) 



(2.9) 



92= T K (r)sin0, p 2 = C 2 (r) sin <j)p r + 

which generalizes (2.1). The inverse change of coordinates leads to the expressions 

C«(r) = -= J==, S K (r) = — ^L=, T K (r) = |q|, 
a/1 + Kq^ a/1 + ftq 2 

J. 91 ■ , 92 , / 92 

cos(p= 1 — r, sm<p=-. — r, tan0= — . 

q q 9i 



(2.10) 
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3 The anisotropic Higgs oscillator 

In order to generalize the integrable anisotropic Hamiltonian on E 2 (2.2) to the sphere and the 
hyperbolic plane we recall that, in terms of geodesic polar coordinates, the Higgs oscillator 
potential on both spaces is given by T 2 (r) [11, 12, 13] (i.e. tan 2 r on S 2 with k = 1 and 
tanh 2 r on H 2 with k = — 1). Furthermore, the superposition of the Higgs oscillator with two 
centrifugal potentials (i.e. the curved SW system) is also known in these variables and defines 
a superintegrable Hamiltonian. Such a system will be our starting point, whose properties 
are summarized as follows. 

Proposition 1. [11, 12, 13] The Higgs oscillator Hamiltonian on S 2 and H 2 with two 
Rosochatius terms (curved SW system) is given in terms of geodesic polar variables by 




At As 

SUr)) ' "~ KV/ ' Sl(r)^ Sf(0"sinV 



U sw = -[p 2 r + 1 + ^ 2 (r) + o2 + „,, , , ■ (3.1) 



The Hamiltonian Hsw Poisson-commutes with the functions 



hi = ~J$l + $ T 2 (r) COS <f> + 

2 T^(r)cos^0 

^02 = J J 2 2 + 5 T 2 (r) sin 2 + (3.2) 
2 T^(r)sm 2 



^r2 + Aitan 2 ^+^ 



where Jij are the generators (2.5). Each set (Hsw, hi, I12) and (%sw> h2, ^12) is formed by 
three functionally independent functions. Moreover, Hsw = hi + h2 + KI12 + k(Ai + A2). 

Essentially, our task now consists in the introduction of an anharmonic fi-term in the 
Hamiltonian (3.1) by keeping the integrability of the system. This process requires to modify 
simultaneously the integrals of motion (3.2). The final result is presented in the following 
proposition, that can be proven through direct computations. 

Proposition 2. The anisotropic Higgs oscillator Hamiltonian on S 2 and H 2 with two 
Rosochatius terms defined in geodesic polar variables as 

f 5 T 2 (r) + n S 2 (r) sin 2 <f> + o2/ ^ , ± + o2 * 2 . 2 , ( 3 - 3 ) 




S 2 K (r)J • ' - ■ S 2 (r)cos 2 ' S 2 (r)sin 2 

Poisson-commutes with the functions 

= \joi + 5 T 2 K (r) cos 2 + -^-^ 



2 ' T 2 (r) cos 2 4> 

T -2,n = J ( Jl 2 + ^12) + $ T 2 (r) sin 2 + S 2 (r) sin 2 <f> + 



r)V UZ 1 ■ ™ Li) ' K\ I "~ T ' KV I T 1 rp2 I \ ■ 2 JL 

2 T^(r) sin <p 



cos 2 ^ sin 2 i 

where hj are the generators (2.5). The pairs (H K ,Ii )Ki ) and (H k ,I2,k) are constituted by two 
functionally independent functions, and H K = + 12,«- 
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(c) (d) 

Figure 2: Some trajectories of the Hamiltonians H K (3.3) and H' K (4.1) on the sphere S 2 with k = 1 and 
without centrifugal terms (Ai = Az = 0). They are plotted in R 3 with ambient coordinates (xo,x) such that 
XQ + xf + X2 = 1. Time runs from t — to t — 8 for the initial data qi = 1, q\ = 1, g2 = —0.5, q2 = 2. (a) the 
superintegrable case of H K with 5 = 1 and S7 = (Higgs oscillator); (b) the integrable H K with 5 = 1 and 
= 0.2; (c) the integrable H K with 5 = 1 and f2 = 3; and (d) the superintegrable W' K with 5=1. 



We stress that for a generic value of Vt the anisotropic Higgs oscillator % K (3.3) is an 
integrable (but, presumably, non superintegrable) O-anharmonic generalization of %sw (3.1) 
since the former is endowed with a single Quadratic integral of motion (either T\ K or ^2,k) 
and its flat k — > limit is not superintegrable. We remark that in this generalization only 
Zi A = Iqi "survives" from the set of three integrals (3.2), meanwhile the two remaining ones 
are mixed within X2 A , and when = we get I^,k = /02 + + K (Ai + A2). Note also that 
the k = Euclidean limit is naturally included in (3.3) and (3.4), and provides (2.2). 

The system defined by proposition 2 can also be written in terms of Beltrami variables by 
applying the transformation (2.10), and it contains only polynomial and rational expresions 

U K = \ (1 + .q 2 ) (p 2 + «(q • p) 2 ) + ^ + il + (1 + .q 2 ) + 

ll,K = \jll + S<ll + h, (3.5) 
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Figure 3: Some trajectories of the Hamiltonians H K (3.3) and T-L' K (4.1) on the hyperboloid H 2 with a = — 1 
and without centrifugal terms (Ai = A2 = 0). They are plotted in R 3 with ambient coordinates (xo,x) such 
that x\ — x\ — x| = 1. Time runs from t = to t = 12 with the initial data qi = 0.1, qi = 0.1, §2 = 0.3, 
q-z = —0.2. (a) the superintegrable case of T-i K with 5=1 and Q = (Higgs oscillator); (6) the integrable T-i K 
with 5 = 1 and f2 = 0.2; (c) the integrable H K with 5 = 1 and Q. = 3; and (d) the superintegrable H' K with 
<5 = 1. 

where now Jy are the functions given by (2.7). In this way, the curved Higgs oscillator 
potential, T re (r), becomes the simple "flat" expression q 2 and the Euclidean Hamiltonian T-L 
is recovered in Cartesian coordinates, in the form (1.1) and (1.2), by setting k = in (3.5). 

Numerical solutions for some particular bounded trajectories of H K (3.3) are explicitly 
illustrated in figures 2 (2D sphere) and 3 (hyperboloid). They are obtained by fixing the 
values of the parameters $7 and 5 (for the sake of simplicity both A« are set equal to zero) 
and by taking the same initial conditions. In this respect, it is worth remarking that: 

• The main dynamical footprint of superintegrable systems consists in the fact that all 
their bounded trajectories are periodic ones. Indeed, when £1 = Ai = A2 = we recover 
the superintegrable Higgs oscillator, whose periodic trajectories are plotted in figures 
2(a) and 3(a). 

• The integrable Hamiltonian (3.3), arising by introducing the anharmonic potential with 
0^0 (we keep Ai = A2 = 0), corresponds to the anisotropic Higgs oscillator. In this 
case the bounded trajectories shown in figures 2(6) and 3(6) turn out to be nonperiodic. 
Therefore, the transition from the Higgs oscillator (O = 0) to the anisotropic one 
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(f2 7^ 0) can be graphically appreciated by comparing figures 2(a) and 2(6) (see also 
3(a) and 3(6)). 

• We also analyse the particular case of H K with Q = 35 and Ai = A2 = 0, which is a 
curved analogue (k / 0) of the superintegrable Euclidean oscillator 1:2 (k = 0) whose 
flat trajectory is the Lissajous 1:2 curve. For the initial conditions here chosen, it is 
found that the curved bounded trajectories are no longer periodic as it is depicted in 
figures 2(c) and 3(c). 

Therefore, the Hamiltonian T-L K (3.3) seems to be superintegrable only in the $7 = case, 
although we stress that we have not rigurously proven the non-superintegrability of T~L K for 
an arbitrary value of Q. Equivalently, we could say that H K presents a superintegrability 
breaking induced by the non-vanishing curvature of the underlying space. 



4 A superintegrable anisotropic oscillator 

As far as the known superintegrable anisotropic system (2.3) is concerned, its curved coun- 
terpart for any value of the curvature k can be written as follows. 

Proposition 3. [11] The Hamiltonian 

l/ 2i V\ \ , . S 2 (r) cos 2 <P 



2\" r Sl(r)) (l- K S 2 (r)cos 2 0) 

+A5 T 2 (r)sin 2 | Ai 

(l- K S 2 (r)cos 2 0) (l-KT 2 (r)sin 2 ^) 2 S 2 (r)cos 2 



is endowed with three integrals of motion given by 

T , 1 , j2 j_ t2 \ 1 r T 2 (r)cos 2 0(l + KT 2 (r)sin 2 0) Ai 

X M- 2 Wn + ^i 2 )+* {l _ KTl{r)sm 2^ + S 2 (r)cos 2 , 

T> 1 T2+AA TE(r)sin 2 



(4.2) 



-2,K 2^02 -T-*" , 2 . 2 a2> 

z (1 — k T^(r) sin 0j 

, T^(r)cos 2 ^sin(/) tan< 

— <A)l<Al2 + ^0 " " " — -9 — 2A 1 



(1 -KT 2 (r)sin 2 (/)) T K (r)cos<£ 

where Jij are the generators (2.5) and %' K = X\ R + X' 2 K . Both sets (H' K ,I[ K , C' K ) and 
(T-L' k ,X' 2 k , JC' K ) are formed by three functionally independent functions. 

Once again, this system can be more easily expressed by introducing Beltrami variables 
through (2.10). In this way we find that 



n K = -{i + ^){ v + «( q • P ) ) + irT ^ fy + (1 - Kq2){1 _ Kql)2 + \i — ^— , 

1 ^2 , „ T 2N , ,9 2 (l + ^2 2 ) , x (1 + ^ 



*M = 2 ( J 01 + ^ 2 ) + 5 (1 - «^)2 + Al — ' (43) 
^ = 2 J 02 + 45 (1 _%2 )2 » £ « = J01 Jl2 + 25 (1 -Zl)2 " 2Al J 
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Figure 4: The potentials U (3.5) andlY' (4.3) (without centrifugal terms Ai = A2 = 0) in Beltrami coordinates 
on the sphere S 2 with k = 1 and on the hyperbolic plane H 2 with k = — 1. (a) the integrable U on S 2 for 
8 — 1 and Q, = 3; (b) the superintegrable U' on S 2 for 5 = 1; (c) the integrable U on H 2 for 5 = 1 and Q — 3; 
and (d) the superintegrable W' on H 2 for 5 = 1. 

where Jij are now the generators (2.7). In this variables, the Euclidean case with k = leads 
just to the equations (1.3) defining %' in Cartesian coordinates. 

The numerical integration of the trajectories for the superintegrable Hamiltonian T~L' K (4.1) 
gives rise to the Lissajous 1:2 curves on the sphere and the hyperboloid which are shown in 
figures 2(d) and 3(d), respectively. They are worthy to be compared with figures 2(c) and 
3(c). Hence, in contrast with the integrable Hamiltonian % K (3.3) with = 35, the system 
H' K can be regarded as the appropriate curved analogue (k ^ 0) of the anisotropic Euclidean 
oscillator 1:2 (k = 0). We recall that in the Euclidean case H' (1.3) is just the particular 
case of % (1.1) with Q = 35 and A2 = 0, but this is no longer true when the curvature is 
turned on. Consequently, we have shown that the Euclidean 1:2 oscillator admits (at least) 
two different integrable curved generalizations, since T~L' K (the superintegrable one) cannot be 
recovered as a particular case of 7i K with £1 = 35 and A2 = 0. 

The difference between these two systems can also be highlighted by representing the 
corresponding potentials. We write H K = T + U (3.5) and H' K = T + W (4.3), with T given 
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by (2.8), and we represent both potentials U (with Q = 35) and U' in Beltrami coordinates in 
figure 4. Notice that the corresponding trajectories for the potentials shown in figures 4(a), 
4(6), 4(c) and 4(d) are drawn in figures 2(c), 2(d), 3(c) and 3(d), respectively. 

On the other hand, the differences between these two curved anisotropic potentials can 
also be enhanced by rewritting them in terms of the ambient coordinates (xo , x) introduced 
in section 2, namely 

c- x2 ^ 9 ^1 -^2 

U = 5- + VLxl + -4 + 



(1 



2 ~ J2 ~ ^2 ' 



W = *TT^+ A& T*Z ff\ 2V2+ ~2> ( 4 - 4 ) 

and recall that xo — > 1 when n — > 0. In this form, these potentials can be more easily compared 
with the proposals for integrable anharmonic oscillators on the spheres and hyperbolic spaces 
that have been previously presented in [18, 19], thus concluding that the anisotropic Higgs 
oscillator potential U defines a new integrable system. In fact, the potentials studied in [18, 
19] are just polynomials when written in ambient coordinates. For instance, the integrable 
potential on the n-dimensional sphere considered in [19] reads 

2 



V = J> 



2 2 




where are arbitrary real constants and Xj are (n + 1) ambient coordinates such that 
Y17=o x f = 1- The systems worked out in [18] have a similar form. Evidently, they are 
quite different from (4.4). 



5 Curved "centrifugal" terms: geometry and dynamics 

It is worth recalling that the SW system on S 2 was interpreted in [20, 21, 22] as a superposition 
of three spherical oscillators, and this result has been extended to higher dimensions and to 
other spaces with constant curvature in [12, 13, 23]. The crucial point is to take into account 
that the Higgs oscillator potential given by T^(r), is a central harmonic oscillator, whose 
center is located at the origin O, and r is the distance between Q and O measured along the 
geodesic I that joins both points (see section 2). Hence if, instead of O, we take a generic 
point Q p such that p is the distance between Q and Q p measured along the geodesic L that 
joins both points, then T 2 (p) can be interpreted as a noncentral harmonic oscillator potential 
with center at Q p . With this ideas in mind, let us consider the potentials U and W written 
in geodesic polar coordinates (3.3) and (4.1) in order to search terms of the type T 2 (p). 

Consider the two base geodesies l\ and li orthogonal at the origin O and the point Q(r, <f>). 
Let Q\ (resp. Q2) be the intersection point of I2 (resp. l±) with its orthogonal geodesic l[ 
(resp. 1' 2 ), x = QQ\ the distance measured along l[ and y = QQ2 the one measured along 
1' 2 as depicted in figure 1. Next we apply the sine theorem [16] on the orthogonal triangles 
OQQi (with external angle <fi) and OQQ2 (with inner angle <f) which yields 

OQQ 1 : S K (x) = S K (r)cos<A; OQQ 2 : S K (y) = S K (r) sin<£. 
Therefore the integrable potential U (3.3) can be expressed as 

U = 6Tl(r) + nSl(y) + + (5.1) 

S K (x) S K (y) 



12 



which provides a unified description of this potential on the three Riemannian spaces: 

• The 5-term is a central harmonic oscillator with center at the origin O. 

• The fi-term is the anharmonic potential which only depends on the geodesic distance y 
to the point Q2, which is not a fixed center since it is moving together with the particle 
located at Q. 

• The Aj-terms (i = 1,2) are two "centrifugal barriers". 
Likewise, the superintegrable potential li' (4.1) can be rewritten as 

U' = ST 2 K (x) + 45T 2 K (y) ( C '( r )C^) + -A-, (5.2) 

\Cl(x)(Cl(r)- K Sl(y)) 2 ) S 2 (*)' 

showing that this potential is not a simple addition of some anharmonic potential (only 
depending on the variable y) to the central harmonic oscillator; its interpretation is as follows: 

• The first term is a noncentral harmonic oscillator with center at the point Qi (which 
is, moreover, not fixed). 

• The second one is not any kind of curved harmonic oscillator, although it could be 
regarded as some kind of "deformation" of a noncentral oscillator with center at Q2- 

• The Ai-term is a "centrifugal barrier". 

It is worth stressing that the Euclidean limit is well defined in both (5.1) and (5.2), namely 



k^O: U = 5r 2 + n y 2 + ^ + ^, U' = Sx 2 + 4Sy 2 + ^ , (5.3) 

x z y z x z 



where now the geodesic distances x and y become Cartesian coordinates (that is, q\ and g 2 
appearing in expressions (1.1) and (1.3)). 

Furthermore, on the sphere S 2 the potential U admits a second interpretation, which 
requires to rewrite the "centrifugal barriers" as noncentral harmonic oscillators. Explicitly, 
let us consider the potential U in the form (5.1), set k = +1 for simplicity and let 0% be the 
fixed point placed along the geodesic k (i = 1,2) which is at a distance | from the origin 
O. Notice that in ambient coordinates (xo,x) we have that O = (1,0,0), 0\ = (0, 1,0) and 
O2 = (0, 0, 1). We introduce the distance ri between Q and Oj measured along such that 
the three points {OiQQi} lie on the same geodesic l[ as it is shown in figure 5(a) in the 
ambient space M 3 . Thus we obtain that 

r\ + x = V2 + y = — , xo = cos r, x\ = sin x = cos n , X2 = sin y = cos T2 ■ 
Then the potential (5.1) on S 2 can be expressed as 

hi = 5 tan 2 r + sin 2 y -\ ^ 1 

sin x sin y 

= Jtan r + — k — r + Ai tan n + A2tan r 2 + Ai + A2- (5.4) 

(l + tan 2 r 2 ) 

Hence this corresponds to the superposition of the central (Higgs) oscillator with center at 
O with an anharmonic potential depending on the distance r 2 to O2 and with two more 
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Figure 5: Geodesic distances involved in the geometric description of the integrable anisotropic Higgs oscil- 
lator with two Rosochatius terms U: (a) on the sphere for the expression (5.4) and (b) on the hyperboloid in 
the form (5.5). 



noncentral oscillators determined by the geodesic distances r« from the centers placed at 
the fixed points Oj. Notice that the flat limit k = can only be performed in the former 
expression of (5.4) providing (5.3) with the two Euclidean centrifugal barriers, since in the 
latter the distances — > oo (that is, the centers Oj go to infinity) so that the terms tan 2 r; L 
become indeterminate. 

We remark that in the hyperbolic plane with k < a similar description in terms of 
oscillators is, in principle, precluded. The analogous points to the previous centers 0- L would 
be beyond infinity as it can be seen in figure 5(6), that is, beyond the "proper" hyperbolic 
space H 2 . In particular, if we set n = — 1 we find that 



xo = coshr, x\ = sinhx, 
and the potential (5.1) on H 2 reads as 

Ai 



x 2 = sinhy, 



U = 5 tanh 2 r + VL sinh 2 y + 



S tanh 2 r + $7 sinh 2 y + 



+ 



A 2 



sinh 2 x sinh 2 y 
Ai A 2 



+ 



tanh x tanh y 



+ Ai + A 2 . 



(5.5) 



Hence the Aj-terms only admit an interpretation as centrifugal barriers with Euclidean limit 
given by (5.3). Nevertheless, an interpretation of the Aj-potentials as noncentral oscillators 
can be established if one considers the complete hyperbolic space. In that case, the centers 
Oi would be located in the "ideal" (exterior) region of H 2 so beyond the "actual" H 2 . 



5.1 Some trajectories 

In order to illustrate the above results and to highlight the role of the curved centrifugal terms, 
we plot in figure 6 some trajectories of the Hamiltonians % K (3.3) and %' K (4.1) on the sphere 
S 2 in ambient coordinates and with a single centrifugal term (Ai > and A 2 = 0). These 
trajectories are numerically integrated by imposing the same initial conditions and values of 
the constants 5 and Q as in the trajectories plotted in figure 2 (which is the same system 
with no centrifugal term). Consequently, the effect of the Ai-potential becomes apparent by 
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(a) 



(&) 




Figure 6: Some trajectories for the Hamiltonians H K (3.3) and H' K (4.1) on S 2 (k = 1) with a single 
Rosochatius potential such that Ai =0.1 and A2 = 0. Time runs from t = to t = 8 with the initial data 
qi = 1, qi = 1, q-2 = —0.5, qy, = 2. (a) the superintegrable case of T-L K with 5 = 1 and Q, = (Higgs oscillator 
with one centrifugal term); (b) the integrable H K with 5=1 and Q = 0.2; (c) the integrable with 5 = 1 
and f2 = 3; and (d) the superintegrable H' K with 5=1. 

comparing figures 2 and 6. In this respect, it has to be taken into account that the Ai-term 
is associated to the ambient coordinate x\ (see (4.4)), with center 0\ = (0,1,0) for the 
noncentral oscillator and geodesic radial distance r\. In this way, we find that: 

• The trajectories are restricted to the "right" semisphere due to the infinite barrier 
arising at x\ = 0. Roughly speaking, this barrier induces a "folding" of the initial 
trajectories in figure 2. 

• The superintegrable Higgs oscillator ($7 = 0) with a centrifugal term provides a kind of 
Lissajous 1:2 curve (figure 6(a)). 

• The integrable anisotropic systems % K with Vt ^ represented in figures 6(6) and 6(c) 
show bounded trajectories which are again nonperiodic. 

• The superintegrable H' K with a centrifugal term/noncentral oscillator gives rise to a 
kind of Lissajous 1:1 curve (figure 6(d)). 

Next, if the second centrifugal term with A2 > is also considered, then similar features 
arise from the second infinite barrier. For the sake of completeness, figure 7 contains the 
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(a) 




(c) 

Figure 7: Trajectories of the Hamiltonian (3.3) on S 2 (k = 1) with two Rosochatius potentials given by 
Ai =0.1 and A2 = 0.05. Initial data are the same as in figures 2 and 6. (a) the superintegrable case of H K 
with 5—1 and Q, = (Higgs oscillator with two centrifugal terms or curved SW system); (6) the integrable 
H K with S — 1 and £1 = 0.2; and (c) the integrable H K with 5=1 and Q = 3. Note that "H^ is not defined 
when both Ai and A2 are non-zero. 

trajectories on the sphere of the Hamiltonian % K (3.3) which exactly correspond to those 
plotted in figure 6 but now with the second Rosochatius potential. We remark that: 

• The trajectories of T~L K are restricted to a quadrant of the sphere; these can be seen as 
once "folded" with respect to the trajectories shown in figure 6 but twice with respect 
to those obtained in figure 2. 

• The superintegrable Higgs oscillator (f2 = 0) with two centrifugal terms (curved SW 
system) gives rise to a Lissajous 1:1 curve (figure 7(a)). 

• As expected once more, the trajectories coming from the integrable anisotropic Hamil- 
tonian ~H K with $7^0 drawn in figures 7(6) and 7(c) show bounded nonperiodic trajec- 
tories. 

Summarizing, we can say that the addition of the appropriate curved centrifugal terms 
preserves the essential dynamical features of all these systems. 



16 



6 Concluding remarks 



In this work we have introduced, ci new integrable Hctmiltonian 1~L^ on S 2 and H 2 , that 
can be interpreted as the superposition of the anisotropic Higgs oscillator with two curved 
Rosochatius potentials. Furthermore, such a system together with the previously known 
curved superintegrable Hamiltonian W K have been algebraically and geometrically studied. 
Numerical solutions for the trajectories corresponding to some particular initial conditions 
for both Hamiltonians have been obtained. In this respect, some remarks are in order. 

• The flat contraction k — > is always a well defined and smooth limit in all the ex- 
pressions presented in the paper. In this way the known systems presented in section 
1 can always be recovered in the Euclidean limit. Hence the two different integrable 
Hamiltonians % K and %' K can be thought as generalizations of the known Euclidean 
systems (1.1) and (1.3) (or (2.2) and (2.3)), and H' K=0 is nothing but a particular 
(superintegrable) case of % k =q. 

• Alternatively, this approach can be seen as a "curvature-deformation" process of the 
Euclidean anharmonic oscillator that by starting from (2.2) and (2.3) arrives at the 
Hamiltonians (3.3) and (4.1), respectively. In this sense, we remark that there is no an 
apparent straightforward prescription in order to perform such a deformation, neither 
in the integrable nor in the superintegrable system, as the explicit expressions for the 
curved Hamiltonians show. 

• It is also worth stressing that when the superintegrable Hamiltonian (4.1) cannot 
be recovered from the integrable one (3.3) by setting 0, = 35 and A2 = (see also the 
potentials (4.4) as well as the integrals Ii ;K (3.4) and X' iK (4.3)). This, in turn, means 
that the Hamiltonians % K and %' K are two completely different systems on the sphere 
and on the hyperbolic plane. This fact indicates that different integrable curved systems 
may collapse under contraction into a single flat one and, conversely, that the integrable 
generalization of a known Euclidean Hamiltonian could split into several generalized 
"curved" possibilites. 

Starting from all these results, several further research directions can be explored. Firstly, 
the search for superintegrable curved generalizations of other commensurate Euclidean oscil- 
lators could be faced. In this respect, we recall that the superintegrable Hamiltonian %' (1.3) 
is the only system in the series of anisotropic oscillators with commensurate frequencies that 
is endowed with quadratic integrals of motion (see [6]). Moreover, the "additional" integral 
with ensures superintegrability is, in fact, a second order reduction of a true six-th order 
integral. More explicitly, if we set Ai = then it is known [6] that 

Qi = (p\ - 25q\ - 2iV2Sq lPl S j 2 (p\ - 85g| + &V2~5q 2 p 2 ) , Qi, Q = \ (Qi + Qi) , 

are sixth order integrals of motion of H', and that the set I[, I' 2 (1-3) and Q is formed by 
three functionally independent functions. Alternatively, the six-th order integral Q can be 
written as 

Q={pj + 25qlf [p\ + 85q 2 2 ) - 166 {26qjq 2 + (qm ~ q2Pi)pif • 

Now, the second order integral C (1.3) arises through the following relation 

,2 _ KL'^X' 2 - Q 
16(5 ■ 
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Remarkably enough, for the curved superintegrable Hamiltonian H' K (4.3) (in Beltrami vari- 
ables), the corresponding six-th order integral Q K can also be found, namely, 



Q . . ( 4l + U ?2 + 2, fi±f) 2 (jg, + ^) - 1M (^ + ^) ^ , 

where J^- are the generators (2.7). Moreover, the relationship with the second order integral 
C' K (4.3) remains in the same "flat" form 

~ 165 " 

Therefore these latter expressions seem to point out that the construction of the S 2 and H 2 
analogues of some other Euclidean commensurate oscillators could be feasible. 

Secondly, another interesting open problem would be the construction of the constant 
curvature analogues of some integrable Henon-Heiles systems (see [24] and references therein) 
that can be written as particular cases of the multiparametric family 

H = \{p\ + vl) + Sqj + (6 + n)q 2 2 + a (qjq 2 + /3(? 2 3 ) + ^, 
z q 1 

in which the only known integrable cases are the Sawada-Kotera (/3 = 1/3,0 = 0), the KdV 
(/3 = 2, fi arbitrary) and the Kaup-Kupershmidt (/3 = 16/3,0 = 155) systems. The results 
here presented would be a very reasonable starting point, that should be completed with 
appropriate curved analogues of the cubic Henon-Heiles term. 

Finally, the explicit solution of the Schrodinger equation associated to Ti K and %' K should 
be addressed, as well as the analysis of the quantum nonlinear dynamics generated by these 
new class of integrable nonlinear quantum models by following, e.g., the approaches presented 
in [25, 26, 27]. Work on all these lines is in progress. 
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